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Abstract 

We propose a Ginzburg-Landau theory for the elastic properties of shape memory polycrystals. A single 
crystal elastic free energy for a system that undergoes a square-to-rectangle transformation is generalized 
to a polycrystal by introducing a crystal orientational field that is determined from a continuum phase 
field model. The coupled system is used to study domain morphology and mechanical properties of shape 
memory alloys in different temperature regimes. 
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I. INTRODUCTION 



Materials known as martensites undergo first-order, diffusionless structural transformations 
from one crystal phase (austenite) to another, usually a twinned, phase (martensite) due to shear 
strains. A subclass of these materials with displacive transformations exhibits the shape memory 
effect^ This effect refers to the existence of a residual strain upon unloading that can be recovered 
on heating to the high temperature austenite phase. In the austenite phase, these materials exhibit 
pseudoelastic behavior where a plateau in the stress-strain curves is observed; however, there is 
no residual strain as the macroscopic deformation is completely recovered when the load is re- 
moved. Interestingly, these unusual mechanical properties of shape memory alloys do not involve 
any plastic effects such as those caused by dislocation motion but are entirely due to the intrinsic 
elastic nonlinearities. As a consequence of these unique properties, shape memory materials have 
many technological applications^. Some examples include NiTi, FePd, AuCd and copper-based 
ternary alloys, e.g. CuAuZn. 

Most commercial applications of shape memory alloys make use of polycrystalline specimens 
and therefore it is important to compare the mechanical response of polycrystals to that of single 
crystals. The problem of finding the effective properties of martensitic polycrystals has been 
studied by analytical methods^*^ and finite element simulations^. However, these methods do not 
account for the complex polycrystal geometry and also do not incorporate the long-range elastic 
interactions between the grains. Continuum simulations that span a range of length scales are good 
candidates to describe these issues. Recently, phase-field micro-elasticity models that employ 
static grains created by the Voronoi construction have been studied^. However, it is important 
to regard the grain orientation as a thermodynamic variable since a polycrystal, in reality, is a 
metastable state that is formed by a grain growth process. In the present work, we attempt to 
account for this metastability of a polycrystal. 

The evolution of grains during grain growth has been studied using the phase-field 
approacbSiSilfl. Although these models correctly describe the grain morphologies and domain 
coarsening, they do not address the issues of elasticity and crystal symmetry. A coupling of these 
models with continuum elasticity models of martensitic transformations provides a framework to 
model mechanical properties of shape memory polycrystals. Here we propose a model in which 
elastic strains are coupled to a phase-field model through an orientation field that is determined 
from a multi-component order parameter describing the crystal orientations. Due to this cou- 
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pling, the strains as well as the grain orientations can change under an external load. In Section 
2 we describe our Ginzburg-Landau model and present two-dimensional (2D) simulations of the 
loading-unloading cycle for martensitic poly crystals in different temperature regimes in Section 3. 
Our main findings are summarized in section 4. 



II. GINZBURG-LANDAU MODEL 

A Landau theory for shape memory materials was first proposed by Falk— . This one dimen- 
sional model captures the salient physics of a number of experimentally observed features of 
martensites. However, the model does not incorporate the elastic long-range interactions that 
are crucial to describe the microstructure of martensites. Barsch and Krumhansl introduced 
a Ginzburg-Landau model to describe the inhomogeneous microstructure of these materialspi^ 
and the model has been extended to simulate martensitic domain structures in two and three 
dimension o ^ 3^ 1 4^ 1 5^ 1 6^ 1 7 n^^t we generalize this theory to describe a 2D square to rectangle marten- 
sitic transformation in a poly crystal with different crystallographic orientations. The theory is 
formulated in terms of a free-energy functional 

F dv grain ~l~ f elastic ~l~ fload\ 5 (1) 

where f grain is the free energy density due to the orientational degrees of freedom of the poly- 
crystal, /elastic represents the elastic free energy and fioad is the free energy contribution due to an 
external applied load. The polycrystalline system is described by a set of Q non-conserved order 
parameters2(7yi, 7^2, •••^Vq)- terms of these order parameters the free energy f grain is given by 



Ua.n = E l^rn' + fry.^ + fry^j + f E E + E f (V^.)^ (2) 

i=l i=l j>i i=l 

For ai, a2 < and a^.a^ > 0, the first two terms in equation (2) describe a potential with Q 
degenerate minima (/yo,0, ...,0), (0,7yo, ....,0) up to (0,0, ....,7yo) (^o > 0), corresponding to Q 
grain orientations. This form differs from that used by Chen and Yang^ as there is an additional 
cubic term in the free energy. The cubic term ensures that r^o ^ and allows us to uniquely 
associate an orientation with each minimum, given by an angle 



(3) 
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Thus, there are Q orientations between and a maximum angle 9^ (as an example, the Q = 3 
case has minima at (ryo, 0, 0) and (0, 770, 0) and (0, 0, ryo) corresponding to 6^ = 0^, 6^ = and 

The gradient energy {K > 0) represents the energy cost of creating a grain boundary. To 
describe elastic effects, we consider the linearized strain tensor in a global reference frame defined 
as Cij = {uij + Uji)/2 [i = 1, 2]j = 1, 2), where ui represents a component of the displacement 
vector and uij is a displacement gradient. We use the symmetry-adapted linear combinations of 
the strain tensor— defined by 61 = {e^x + ^yy)lV2, 62 = {e^x — ^yy)/V2 and 63 = e^y. Here ei 
represents the bulk (dilatation) strain, 62 the deviatoric (rectangular) strain and 63 the shear strain. 
To generalize these definitions to the case of a poly crystal described by an orientational field 9{fj), 
the strain tensor in a rotated frame is calculated using R{9{ff)) e R^(9{ff)), where R{9{ff)) is a 
rotation matrix and 9{ff) is determined from the minima of Fg^ain using equation (3). Under this 
rotation, the elastic free energy describing a square to rectangle transition is given by 

J elastic = + —63 + flocal[e2) + —[^62) , (4) 

where ei = ei, 62 = 62 cos [2(9 (ry)] + ^263 sin[2(9(77)] and 63 = -(1/a/2) 62 sin [2(9(7/)] + 
63 cos[26^(f/)]. The nonlinear part of the elastic free energy is 

fiocaiie2) = ^e,' + je,^ + ^e,^ (5) 
z 4 o 

At the level of the unit cell, the transformation free energy fiocai describes a square to rectangle 
transition where the austenite phase (square) has 62 = and the martensite phase is described by 
62 = ±eo corresponding to the two rectangular variants. Here Ai = Cu + C12, A2 = Cu — 
C12 and As = 4C44, where Cu, Cu and C44 are the elastic constants for a crystal with square 
symmetry. The quantities a and /? are the higher order nonlinear elastic constants and K2 is 
the appropriate deviatoric strain gradient coefficient that determines the energy cost of creating 
a domain wall (twin boundary) between the two rectangular variants. In principle, K2 can be 
determined experimentally from phonon dispersion curved. In this work, we are interested in 
simulating a uniaxial loading experiment. If we choose the x axis to be the loading axis, then the 
free energy contribution due to the external load a is given by 

fload —(^^xx ~^ 



a 

71 



ei + 62 cos[26'(?7)] - V2e3 sm[2e{ff)] 



(6) 



The strains ei , 62 and €3 are not independent but are related by the elastic compatibility relationi^ 



V% - (- 



)e2 - VStt^cs = 0. 



' dx"^ dy"^ dxdy 
In terms of the strain variables ei, 62 and 63, this relation becomes 



(7) 
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62 cos[2^(?7)] - ^263 sin[2^(?])] 



dxdy 



V2 



0. 



(8) 



We eliminate the strain ei from f elastic and fioad using equation (8) so that the total elastic energy 
is in terms of 62 and 63 only. This method has been used previously in the context of martensitic 
transformations in single crystal^. We first introduce e[ = ei — {a/Al^/2) so that /e// = 

f elastic + fload. whcrC 



,2 , ^3 2 



feff = -^ei + ^^3 + f local {e2) ^ 



62 cos[29{f])] - V2essm[29{f])] 



(9) 



To eliminate e[ from equation (9), we use the Fourier representation of equation (8) to obtain 



kx H~ k' 



^/8kx ky 
kx H" ky 



(10) 



Vy IVx I IVy 

where T2(k) represents the Fourier transform of 62 cos[26^(f/)] — V2e^ ^m[29{ff)\ and T^{k) is the 
Fourier transform of e2(sin[2^(rf)]/\/2) + 63 cos[2^(7/)]. Thus, the effective free energy can be 
written as 



eff 
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dk 



kx H" ky 



2\ 2 
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y/8kxky(^kx ky ) 
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T,{k)T2{-k) + r3{-k)r2{k) 



df 



^63' + fiocai{e2) - (^62 co8[29{f])] - v^gg sin[2e(7y)] 



(11) 



The long-range part of the free energy is always orientation dependent due to the elastic com- 
patibility induced anisotropy in the kernels. This term ensures that elastic compatibility will be 
satisfied within the grains as well as at the grain boundaries. We emphasize that the long-range 
interaction excludes the |fc| = mode since compatibility is trivially satisfied for this case. The 
|fc| = mode, which refers to the homogeneous state, is accounted for by the local terms in F^ff. 
The total free energy of the system is F = Fgrain + Feff^ and we assume relaxational dynamics 
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for 62 and 63, that is, 



de2 




'5F' 
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_Se2_ 






'SF' 


dt 


-73 


_5e3_ 



(12) 

where 72 and 73 are the appropriate kinetic coefficients for the deviatoric and shear strains. Simi- 
larly, the dynamics of the grains is defined by Q equations 

5F 



dt 



5m 



(13) 



where 7^ is a kinetic coefficient and i = 1, Q, corresponding to Q grain orientations. 



III. SIMULATIONS OF TEXTURE AND STRAIN EVOLUTION 

We simulate the mechanical properties of shape memory materials using the model described 
in Section 2. We choose FePd parameter^ for which Ai = 140 GPa, As = 280 GPa, a = 
— 1.7 X 10^ GPa and /? = 3 x 10^ GPa. The temperature dependent elastic constant A2 undergoes 
a softening and hence controls the square to rectangle transformation. We study four different 
cases corresponding to A2 = —2, 1, 2, 3 GPa. Muto et al.— have measured the elastic constants 
of FePd as a function of temperature. The lowest temperature measurement they reported was 
at 290 K corresponding to ^2 ^ 10 GPa. Thus the values of A2 we have chosen correspond 
to temperatures lower than 290 K and are in the vicinity of the transition temperature of 265 K. 
Figure 1 shows the profiles for the free energy fiocai for these values. For the parameters in f grain 
we choose (for illustrative purposes) ai = — 10 GPa, a2 = — 10 GPa, = 10 GPa, a4 = 20 GPa, 
Q = 5 and 9^ = 30^. Here, we also need to specify the grain boundary energy coefficient K and 
the strain gradient coefficient K2. For FePd, the strain gradient coefficient^'^ has been measured to 
be K2/a^^ = 25 GPa, where ao is the lattice spacing of the crystal. The grain boundary energy 
coefficient is chosen as Kja^ = 10^ GPa. The lengths are scaled by f = (lOOao)C- For a 
homogeneous single crystal, using these parameter values, the free energy in equation (2) has 
5 degenerate minima defined by 6^0 (^) = 0^, 7.5'', 15"", 22.5^, 30^. We should point out that we 
can consider a more continuous orientation distribution of the polycrystal with a large number of 
states (large value of Q) in the full range (0"" — 45^). However, for the sake of clarity of the domain 
patterns, we restrict ourselves to the range described above. Equations (12) and (13) are solved 
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numerically to simulate the domain structures and mechanical properties in different regimes. For 
simplicity, we assume 7ry = 72 = 73 = 7 and use rescaled time defined by t* = t(10^°7). 

An initial poly crystalline configuration is first generated by solving equations (12) and (13) for 
a 12800ao x 12800ao system with periodic boundary conditions, starting from random initial condi- 
tions. A grain growth process is simulated with a = in the austenite phase so that all components 
of the strain tensor vanish. Grains with orientations ^o(^) = 0^, 7.5^, 15^, 22.5^, 30'' form and start 
coarsening. We arrest the system in a given polycrystalline configuration by abruptly changing 
the value of the parameter ai from —10 GPa to —160 GPa (the parameter A2 is also changed so 
that the system is in the desired martensitic phase). This sudden decrease in ai increases the free 
energy barriers between the crystalline states and the growth stops. We first consider the case 
A2 = —2 GPa. Figure 1 shows that for a homogeneous system, the transformation free energy 
f local {^2) for A2 = —2 GPa has a local maximum at 62 = and two degenerate global minima. 
In the absence of applied stress (a = 0), the arrested polycrystal evolves into a domain pattern of 
the variants of the martensitic phase (there is no austenite present since 62 = is unstable). The 
strains in each grain as well as the orientation of the martensitic domain walls (i.e. twin bound- 
aries) are determined by the orientation of the grain. This behavior is clear from Figure 2(a) that 
shows the distribution of the strain 62 (r) (deviatoric strain relative to the global frame of reference) 
and Figure 2(b) that shows the local orientation 9{f]{r)). The domain walls are oriented at angles 
6{f]) + 7r/4 or 6^(f/) — 7r/4. We point out that the average strains for this configuration are very 
small and correspond to a system with no macroscopic deformation. 

To simulate mechanical loading, an external tensile stress a is applied quasi- statically, i.e., 
starting from the unstressed configuration of Figures 2(a) and 2(b), the applied stress a is increased 
in steps of 5.13 MPa, after allowing the configurations to relax for t* = 25 time steps after each 
increment. The loading is continued till a maximum stress of a = 200 MPa is reached. Thereafter, 
the system is unloaded by decreasing a to zero at the same rate at which it was loaded. Figures 2(c) 
and 2(d) relate to a stress level of a = 46.15 MPa during the loading process. The favored variants 
(red domains in the left panel) have started to grow at the expense of the unfavored variants (blue 
domains in the left panel). The orientation distribution 9{f]{f)) in Figure 2(d) has not changed 
much. As the stress level is increased further, the favored variants grow. Even at the maximum 
stress of 200 MPa, some unfavored variants persist, as is clear from Figure 2(e) (in fact, further 
application of stress does not remove such structures as the unfavored variants are in a "locked" 
state due to intergranular constraints). We note that the grains with large misorientation with the 
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loading direction rotate, as is clear by comparing Figure 2(f) with Figure 2(b) and Figure 2(d), 
where we can observe that dark red colored grains have turned orange, indicating rotation of those 
grains. Grains with lower misorientation do not undergo significant rotation. This rotation is due 
to the tendency of the system to maximize the transformation strain in the direction of loading so 
that the total free energy is minimized. Within the grains that rotate, sub-grain bands with slightly 
higher values of the orientation 9{f]{f)) are present. These bands correspond to the unfavored 
strain variants that still survive. Figure 2(g) and 2(h) show the situation after unloading to a = 0. 
Upon removing the load, a domain structure is nucleated again due to the local strain gradients at 
the grain boundaries and the surviving unfavored variants in the loaded polycrystal configuration 
in Figure 2(e). This domain structure is not the same as that prior to loading (Figure 2(a)) and 
thus there is an underlying hysteresis. The unloaded configuration has non-zero average strain. 
This average strain is recovered by heating to the austenite phase, as per the shape memory effect. 
Figure 2(h) shows that the orientation distribution reverts to its preloading state as the grains rotate 
back when the load is removed. To show the rotation of grains more clearly, we plot in Figure 3 the 
change in orientation /\9 = 9 {a) — 9 {a = 0) for a = 200 MPa. It is clear by comparing this figure 
with Figure 2(b) that only grains with large misorientations with the loading axis {9 = 0^) rotate 
in order to decrease the misorientation. In some of these rotating grains, we observe bands where 
the misorientation has increased. As discussed earlier, these bands correspond to regions where 
the unfavored variants do not disappear even at high stresses. This behavior is a consequence of 
the intergranular constraints. 

We compare the above mechanical behavior of the polycrystal to the corresponding single crys- 
tal. A single crystal simulation is set up with exactly the same elastic free energy parameters as 
the polycrystal case but the orientation is fixed at 6^ = 0"". Loading conditions are also identi- 
cal to the polycrystal case. Figure 4 shows the evolution of the martensitic variants during the 
loading-unloading cycle for the single crystal. Figure 4(a), shows the simulated microstructure 
(12800ao X 12800ao) prior to loading, with domain wall orientations at 7r/4 or — 7r/4 everywhere. 
In contrast, the domain wall orientations change from grain to grain in the polycrystal case depicted 
in Figure 2(a). Figure 4(b) shows the domain patterns for a stress a = 46.15 MPa. We see that 
the favored variants (red domains) grow at the expense of the unfavored variants (blue domains). 
This growth continues, as shown in Figure 4(c) (a = 56.41 MPa) and Figure 4(d) (a = 71.79 
MPa). Finally, in Figure 4(e) (a = 200 MPa) all unfavored variants disappear and we obtain a 
single domain, in contrast to the polycrystal case at the same stress level (Figure 2(e)) where unfa- 
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vored variants persist. Upon unloading, the single crystal remains in a single domain state (Figure 
4(f)). This is because there are no inhomogeneities or thermal noise in the simulations to cause 
renucleation of the domain structure and so the system remains in the positive strain minimum 
homogeneously. The unloaded poly crystal, however, reverts to the domain pattern (Figure 2(g)). 
This also shows that, unlike the single crystal case, the mechanical behavior of the polycrystal is 
not necessarily governed by fiocai{^2) but by a more complex inhomogeneous free energy land- 
scape. A consequence of this is that the residual strain for a single crystal^i will be larger than that 
for a polycrystal^. The stress-strain curves corresponding to Figure 2 and Figure 4 are shown in 
Figure 5. The residual strain for the polycrystal (~ 0.7%) is smaller than that for the single crystal 
(^ 1.8%) due to an effective averaging over different orientations and nucleation of domains at 
grain boundaries upon unloading. Also, the change in the stress-strain curve for the polycrystal is 
not abrupt because the response of the polycrystal is averaged over all grain orientations. 

We now discuss cases where the material exists in the austenite phase prior to loading. As 
can be seen in Figure \, A2 = 1,2 and 3 GPa correspond to this situation. For A2 = I GPa, 
the austenite (62 = 0) is a metastable local minimum and there are two stable martensitic minima. 
Figure 6(a) and 6(b) show the situation for a = 0. Since there are no nucleation mechanisms in the 
present simulations, the system remains in an austenite phase (62 = 0) in all the grains. Loading 
is simulated in exactly the same manner as for the ^2 = — 2 case. In Figure 6(c) and Figure 
6(d) we observe that a stress induced transformation has heterogeneously occured in all the grains 
(a = 30.76 MPa) (although no local stresses have been introduced in the model, the heterogeneous 
nucleation occurs due to the strain gradients at the grain boundaries which act as embryos for 
the transformation). The transformed regions are represented by yellow/red shades whereas the 
untransformed austenite regions are represented by light green shades. Some grains also show the 
presence of the unfavored variant (62 < 0) represented by dark green and blue shades. Figure 6(d) 
shows that at this stress level, the texture is the same as for the unloaded configuration. At the 
maximum stress (Figure 6(e)) a = 200 MPa, most of the system has transformed although some 
austenite and the unfavored variant domains remain. As in the earlier case, grains with higher 
misorientation have also rotated (Figure 6(f)). Upon unloading, for this case also some of the 
unfavored domains reappear, as can be seen in Figure 6(g). Figure 6(h) shows that the texture 
returns to its preloading state. 

Next, we consider the case ^42 = 2 GPa for which 62 = is the global minimum and there 
are two metastable martensitic minima. For this case, the "arrested" polycrystal exists in the 
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austenite phase before loading. This situation is depicted in Figure 7(a) where the strain 62 is 
close to zero everywhere. The corresponding orientation distribution is shown in Figure 7(b). 
The system is loaded at an identical loading rate as in the earlier cases. Figure 7(c) shows the 
spatial distribution of 62 at a load of a = 30.76 MPa. A stress induced martensitic transformation 
from the square phase to the rectangular phase takes place heterogeneously, as can be seen from 
the transformed regions (yellow/red shades) that are embedded in a matrix of the untransformed 
austenite (green/blue shades). Figure 7(d) shows that at this stress level, the orientation distribution 
does not change appreciably. Upon loading further to a = 200 MPa, the transformed phase grows 
although the austenite phase is locally retained in some regions (Figure 7(e)). Similar to the earlier 
cases, transformation accommodating grain rotations are also observed for this case, as seen in 
Figure 7(f). When the load is removed, the domain structure shown in Figure 7(g) is observed. 
Here, most of the system reverts to an austenite phase but some martensitic domains (as yellow 
streakings) remain giving rise to a very small residual strain. Figure 7(h) shows the corresponding 
orientation distribution that returns to its preloading state. 

Finally, we discuss the case ^2 = 3 GPa which delineates a pseudoelastic behavior. The 
appropriate free energy shown in Figure 1 has only one minimum, corresponding to 62 = 0. For 
this case as well, the preloading state is austenite and the strain 62 is very small everywhere, as 
seen in Figure 8(a). The corresponding orientation distribution is shown in Figure 8(b). Nucleation 
of martensitic variants in the austenite matrix can be observed in Figure 8(c), corresponding to 
(J = 30.76. Figure 8(d) shows that the orientation distribution associated with this stress does 
not change much. As the stress is further increased to a = 200 MPa, most of the system is 
transformed to (stress induced) martensite, as shown in Figure 8(e). Rotations of grains with large 
misorientation with the loading axis are observed for this case also (Figure 8(f)). In contrast to the 
^42 = —2 GPa and ^2 = 2 GPa cases, this system reverts to a homogeneous austenite phase upon 
unloading (Figure 8(g)), i.e. there is no residual strain. Thus, this case shows pseudoelasticity. 
The grain rotations are also recovered upon unloading, as seen in Figure 8(h). 

We have also simulated the single crystal loading-unloading cycle corresponding to Figures 
6, 7 and 8. Unlike the ^2 = — 2 GPa case discussed in Figure 3, the system always remains in 
a single domain state during the loading and unloading. This is because there are no nucleation 
mechanisms in the single crystal cases to create twin boundaries. A stress induced martensitic 
transformation occurs homogeneously at a critical load, for all these cases. Upon unloading, there 
is no renucleation of the twinned state. 
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The stress-strain curves corresponding to Figures 6, 7 and 8 are depicted in Figure 9. The 
corresponding single crystal curves are also shown. It is clear that the residual strain decreases as 
the elastic constant A2 (or the temperature) increases. It is also seen that the residual strains for 
the polycrystal simulations corresponding to = 1 GPa and A2 = 2 GPa are smaller than their 
single crystal counterparts due to averaging over different orientations and nucleation mechanisms 
that exist in the polycrystal. However, for the pseudoelastic case (A2 = 3 GPa), the residual strains 
are zero for both the single- and polycrystal cases. 

IV. SUMMARY AND DISCUSSION 

We have proposed a framework to study the mechanical properties of shape memory polycrys- 
tals. We have coupled the elastic free energy for a square to rectangle transition to a phase-field 
model describing crystal orientations. This approach can be readily extended to any crystal sym- 
metry and does not require any a priori assumption of grain shapes or microstructure. The mi- 
crostructure is governed by the crystal symmetries encoded in the appropriate elastic free energy 
functional. The long-range elastic interaction between the grains is also incorporated. An impor- 
tant feature of the present work is the coupling between the grain orientation and elasticity so that 
the metastability of the polycrystal may be accounted for within the same framework. 

We studied mechanical properties of shape memory polycrystals and single crystals in different 
temperature regimes. There are significant differences between the mechanical response of single- 
and polycrystals. Since the mechanical properties of the polycrystal are an average of individual 
grains, the stress-strain curves are smoother compared to those of the single crystals. The inhomo- 
geneities in the polycrystal ensure that domain walls influence the mechanical behavior throughout 
the loading-unloading process. In the temperature regimes with nonzero residual strain, the un- 
loaded polycrystals not only have reduced strain but also show domain microstructure. In contrast, 
the simulated defect free single crystals exhibit no such patterns after unloading and have much 
higher residual strains. Our findings are consistent with the fact that in general, polycrystals have 
poor shape memory properties in comparison to single crystals^iifii^l. We emphasize that our find- 
ings apply to materials where the symmetry change due to the transformation is not particularly 
large and the number of low symmetry variants is relatively small (for example materials that show 
a cubic to tetragonal transformation). Materials such as NiTi, which undergo a cubic to monoclinic 
transformation have a relatively larger symmetry change and larger number of low symmetry vari- 
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ants. For NiTi, the differences between shape memory behavior of polycrystals and single crystals 
is not significant. Thus, the issues of how well strain is accommodated and how favored are the 
variants can play a role in determining the shape memory properties. Clearly, further investigation 
of this aspect is required. 

The simulations also predict that grains with a large misorientation with the loading direction 
can undergo reversible rotations, even in the absence of plastic flow. The magnitude of the rotation 
is very small in the linear elastic regime and significant rotation occurs only after the transforma- 
tion begins. It is well known that grains can rotate in the plastic regime to accommodate crystal- 
lographic slip2^. Although there are no plastic effects within the present model, there are elastic 
nonlinearities that couple to the orientation degrees of freedom. Measurements of texture evolu- 
tion upon loading in NiTi shape memory alloys have recently been reported^^. However, from 
these experiments, it is not possible to clarify if the rotational mechanism we observed exists. Our 
simulations suggest that for grains with a high misorientation with the loading axis , in principle, 
rotation as well as detwinning can occur simultaneously for shape memory polycrystals. There 
appears to be an underlying tendency in the poly crystal to decrease the effect of misorientations. 
Based on experiments and finite element simulations, a similar tendency where favorably oriented 
grains induce the transformation in unfavorably oriented grains, thereby effectively reducing the 
effect of grain misorientations, has been discussed by Gall et al^. Since the grain orientations 
are not fixed in our simulations, this reduction in the effect of misorientations can be achieved by 
evolution of the transformation strains as well as the grain orientations. Further experiments are 
needed to clarify to what extent such a mechanism indeed exists in real shape memory polycrys- 
tals. The magnitude of this grain rotation depends on the parameters of the poly crystal phase field 
model, or on the energy barrier between the grains. These parameters have been chosen arbitrar- 
ily for illustrative purposes in the present simulations. An estimation of these parameters from 
experiment or atomistic simulations will be useful for a better comparison with experiments. 
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FIG. 1: Local free energy density vs. 62 as a function of temperature. The relatively high temperature 
(black curve) has A2 = 3 GPa, and successive curves with decreasing temperatures are for ^2 = 2, 1 and 
—2 GPa, respectively. 
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FIG. 2: Spatial distribution of 62 (r) (deviatoric strain in a global frame) for a poly crystal with A2 = —2 
GPa (snapshots (a),(c),(e) and (g)) and 9{f) in degrees (snapshots (b),(d),(f) and (h)). The corresponding 
stress levels are a = (before loading, (a) and (b)), a = 46.15 MPa ((c) and (d)), a = 200 MPa ((e) and 
(f)) and cr = (after unloading, (g) and (h)). 16 



FIG. 3: Spatial distribution of AO ^ 9(^a) — 9{a ^ 0) in degrees for a = 200 MPa (corresponding 
Figure 2(f)). 
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FIG. 4: Spatial distribution of 62 (r) (deviatoric strain in a global frame) for a single crystal {0{f]) = 0^) 
with A2 = —2 GPa. The corresponding stress levels are a = (before loading) (a), a = 46.15 MPa (b), 
a = 56.41 MPa (c), a = 71.79 MPa (d), a = 200 MPa (e) and a = (after unloading) (f). 
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UNIAXIAL STRAIN 



FIG. 5: Variation of applied load a with average uniaxial strain {cxx) for A2 
(red curve) and a single crystal (blue curve). 



19 






FIG. 6: Spatial distribution of 62 (r) (deviatoric strain in a global frame) for a poly crystal with = 1 GPa 
(snapshots (a),(c),(e) and (g)) and 9{f) (snapshots (b),(d),(f) and (h)). The corresponding stress levels are 
(7 = (before loading, (a) and (b)), a = 30.76 MPa ((c) and (d)), a = 200 MPa ((e) and (f)) and a = 
(after unloading, (g) and (h)). 20 





FIG. 7: Spatial distribution of 62 (r) (deviatoric strain in a global frame) for a poly crystal with A2 = 2 GPa 
(snapshots (a),(c),(e) and (g)) and 9{f) in degrees (snapshots (b),(d),(f) and (h)). The corresponding stress 
levels are a = (before loading, (a) and (b)), a = 30.76 MPa ((c) and (d)), a = 200 MPa ((e) and (f)) and 
a = (after unloading, (g) and (h)). 21 




FIG. 8: Spatial distribution of 62 (r) (deviatoric strain in a global frame) for a poly crystal with = 3 GPa 
(snapshots (a),(c),(e) and (g)) and 9{f) in degrees (snapshots (b),(d),(f) and (h)). The corresponding stress 
levels are a = (before loading, (a) and (b)), a = 30.76 MPa ((c) and (d)), a = 200 MPa ((e) and (f)) and 
a = (after unloading, (g) and (h)). 22 



(/) 100 

_l 
< 




FIG. 9: Variation of applied load a vs. average uniaxial strain (cxx) as a function of (increasing) tempera- 
ture. From top to bottom, A2 = 1, 2, 3 GPa, respectively. Red curves correspond to a poly crystal and blue 
curves to a single crystal. 
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